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Abstract 

We present results from a measurement of double diffraction dissociation in pp collisions at 
the Fermilab Tevatron collider. The production cross section for events with a central pseudora- 
pidity gap of width Ar;" > 3 (overlapping ?; = 0) is found to be 4.43 ± 0.02(stat)±1.18(syst) mb 
[3.42 ± 0.01(stat)±1.09(syst) mb] at ^/s — 1800 [630] GeV. Our results are compared with previous 
measurements and with predictions based on Regge theory and factorization. 

PACS number(s): 13.85.Ni 

Double diffraction (DD) dissociation is the process in which two colliding hadrons dissociate into clusters of particles 
producing events with a central pseudorapidity |1| gap (region of pseudorapidity devoid of particles), as shown in 
Fig. ??. This process is similar to single diffraction (SD) dissociation, in which one of the incident hadrons dissociates 
while the other escapes as a leading (highest momentum) particle. Events with pseudorapidity gaps are presumed to 
be due to the exchange across the gap of a Pomeron [EJ, which in QCD is a color singlet state with vacuum quantum 
numbers. 

Previous measurements of DD have been performed only over limited pseudorapidity regions for pp collisions 
at ^/s — 200 and 900 GeV [g|, for exclusive and semi-inclusive dissociation channels at lower energies [^|j§|, e.g 
PP ~^ (p7r+7r^)(p7r+7r^) or pp -^ {p7r^7r^)+X, and for ^p interactions at the DESY ep collider HERA |^. The present 
measurement, based on a study of central rapidity gaps in minimum bias events from pp collisions at ^/s = 1800 and 
630 GeV collected by the Collider Detector at Fermilab (CDF), covers a wide ry range, allowing comparisons with 
theoretical predictions on both ry-dependence and normalization. 

To facilitate our discussion, we begin by defining the relevant variables [Q]. We use s and t for the square of the 
c.m.s. energy and 4-momentum transfer between the two incident hadrons, ^ for the fractional momentum loss of 
the leading hadron in SD, and 77 for pseudorapidity. For pp double diffraction dissociation into masses Mi and M2, 
we define the nominal pseudorapidity gap as At] = In . .^ ° ^ , where So = 1 GeV^; on average, the nominal gap is 

approximately equal to the true rapidity gap in an event. A variable defined as s' = MIM2/ Sq can be thought of 
as the generalization of s' = A/^ for SD, since in both cases In — represents the pseudorapidity region accessible 
to the dissociation products of the diffractive sub-system(s). For pp SD with M2 = nip sa 1 GeV, s' = A/f and 



-Ar; 



s'/s. 



Diffraction has traditionally been treated theoretically in the framework of Regge phenomenology [g[. At large Ar], 
where Pomeron exchange is dominant |l7|, the SD cross section is given by the triple-Pomeron amplitude. 
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where a{t) is the Pomeron trajectory, (3{t) the coupling of the Pomeron to the (anti)proton, and k = g{t)/f3{0) 
the ratio of the triple-Pomeron to the Pomeron-proton couplings; we use a{t) — a{0) + a't = 1.104 -I- 0.25i [pj, 
/3(0) = 4.1 mbi/2 g^ ai^d ^(^) ^ o.69 mh^/^ (^ k = 0.17) ||l. The second factor of Eq. has the form of the 
Pomeron-proton total cross section at the sub-energy VV, while the first factor can be thought of as a rapidity 
gap probability |l^. Measurements on SD have shown that Eq. ^ which is based on Regge factorization, correctly 
predicts the A?7 dependence for Ar/ > 3, but fails to predict the energy dependence of the overall normalization, 
which at -^s = 1800 GeV is found to be suppressed by an order of magnitude |l^,|l^. It is generally believed that 
this breakdown of factorization is imposed by unitarity constraints p3| . Phenomenologically, it has been shown that 
normalizing the integral of the gap probability (first factor in Eq. 1) over all phase space to unity yields the correct 
energy dependence |P,|l2|. 

Using factorization, the DD differential cross section may be expressed in terms of the SD and elastic scattering 
cross sections as [0 
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where boo = 2a' In (sSo/M^M|). Changing variables from Mi and M2 to Ar] and r]c = In ^, where r]c is the center 
of the rapidity gap, yields (setting /3i = /32 => /3) 
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This expression is strikingly similar to Eq. |l], except that, since the gap is now not adjacent to a leading (anti)proton, 
Tjc is treated as an independent variable. The question that arises naturally is whether Eq. ^ correctly predicts the 
differential DD cross section apart from an overall normalization factor, as is the case with Eq. n^for SD. The answer 
to this question, and the suppression in overall normalization relative to that observed in SD, could provide a crucial 
check on models proposed to account for the factorization breakdown observed in SD. 

The components of CDF [|^ relevant to this study are the central tracking chamber (CTC), the calorimeters, and 
two scintillation beam-beam counter (BBC) arrays. The CTC tracking efficiency varies from ^ 60% for pT = 300 
MeV to over 95% for px > 400 MeV within \ri\ < 1.2, and falls monotonically beyond \ri\ = 1.2 approaching zero at 
I77I ~ 1.8. The calorimeters have projective tower geometry and cover the regions \ri\ < 1.1 (central), 1.1 < \ri\ < 2.4 
(plug), and 2.2 < \ri\ < 4.2 (forward). The Arj x A(j> tower dimensions are 0.1 x 15° for the central and 0.1 x 5° for 
the plug and forward calorimeters. The BBC arrays cover the region 3.2 < |?7| < 5.9. 

Events collected by triggering on a BBC coincidence between the proton and antiproton sides of the detector 
comprise the CDF minimum-bias (MB) data sample. Two MB data sets are used in this analysis, consisting of 
1.0 X 10^ ( 1.6 x 10^) events at ^/s = 1800 [630] GeV obtained at average instantaneous luminosities of 2.5 x 10'^° 
(9.6 X 10^®) cm~^sec~"'^. At these luminosities, the fraction of overlap events due to multiple interactions is estimated 
to be 20.7 (6.5)%. To reject overlap events, we accept only events with no more than one reconstructed vertex within 
±60 cm from the center of the detector. 

The method we use to search for a DD signal is based on the approximately flat dependence of the event rate on 
Ar] expected for DD events, as seen by setting a{t) — 1.104 + 0.25t in Eq. ||, compared to the exponential dependence 
expected for non-diffractive (ND) events where rapidity gaps are due to random multiplicity fluctuations. Thus, in 
a plot of event rate versus Ary, the DD signal will appear as the flattening at large A77 of an exponentially falling 
distribution. For practical considerations, our analysis is based on experimental gwps deflned as Arj^^ = rjmax —rjmin, 
where {rjmin) Vmax is the rj of the "particle" closest to 77 = in the (anti)proton direction (see Fig. 1). A "particle" 
is a reconstructed track in the CTC, a calorimeter tower with energy above a given threshold, or a BBC hit. The 
(uncorrected) tower energy thresholds used, chosen to lie comfortably above noise level, are Et = 0.2 GeV for the 
central and plug and E = \ GeV for the forward calorimeters. At the calorimeter interfaces near \ri\ ~ 0, 1.1 and 
~ 2.4, where the noise level is higher, JTyj-dependent thresholds are used. The DD signal is extracted by fitting the 
measured Arj^^^ distribution with expectations based on a Monte Carlo (MC) simulation incorporating SD, DD and 
ND contributions. The same thresholds are used in the MC simulations after dividing the generated particle energy 
by an //-dependent energy calibration coefficient representing the ratio of true to measured (uncorrected) calorimeter 
energy [^. For charged-particle tracks, the MC generation is followed by a detector simulation. 

Figure ?? shows lego histograms of events versus rjmax and —rjmin for data and for Monte Carlo generated ND, SD 
and DD events at ^/s — 1800 GeV. A uniform 77-distribution was assumed for particles within a calorimeter tower. 
The observed structure in the distributions along ?7maa:(mm) is caused by the variation of the tower energy threshold 
with \ri\. The bins at \r]max{min)\ = 3.3 contain all events within the BBC range of 3.2 < \rimax{min)\ < 5.9. 

The diffractive Monte Carlo program is a modifled version of that used in Ref . Mm , incorporating the differential 
cross sections of Eqs. |l| and ||. Non-diffractive interactions are simulated using PYTHIA ||l^. The data distribution 
in Fig. ?? has a larger fraction of events at large \r]max(m.in.)\ than either the ND or the SD Monte Carlo generated 
distributions. From the previously measured SD cross section |l|] and the MC determined fraction of SD events 
triggering both BBC arrays, the fraction of SD events in our 1800 [630] GeV data sample is estimated to be 2.7% 
[2.4%]. A combination of 97.3% ND and 2.7% SD generated events cannot account for the data at large \rimax(min)\ 
in Fig. ??. The simulated DD distribution is approximately flat in \rimax(min)\ and describes the data well when 
combined with the ND and SD distributions, as shown below. 

Figure ?? presents the number of events as a function of Arj^ for the 1800 GeV data (points) and for a flt to the 
data using a mixture of MC generated DD and "non-DD" (ND plus SD) contributions (solid histogram). The dashed 
histogram shows the non-DD contribution. The agreement between data and MC indicates that, as in SD, the shape 
of the differential DD cross section is correctly described by Regge theory and factorization. 

At yi = 1800 [630] GeV, the fraction of events with Ar/^^p > 3 (gap fraction) is (9.08±0.03)% [(11.43±0.03)%], in 
which the fraction of background non-DD events, estimated using the MC simulation, is (23.6 ±0.6)% [(29.7 ±0.6)%]. 
After background subtraction, the DD gap fraction becomes (6.94 ± 0.06)% [(8.03 ± 0.08)%]. The quoted errors 
are statistical. The amount of ND background in the region Ar]^ > 3 depends on the tower energy calibration 



coefficients and thereby on the calorimeter tower energy thresholds used in the MC. Increasing these thresholds has 
the effect of decreasing the multiplicity in the MC generated events, resulting in larger rapidity gaps and hence 
larger ND backgrounds in the region of lS.rfl . The systematic uncertainty in the background is estimated by raising 
(lowering) the tower thresholds in the MC by a factor of 1.25 and refitting the data. This change in thresholds 
increases (decreases) the background by a factor of 1.54 (0.52) [1.56 (0.56)]. 

The vertex cut employed to reject events due to multiple interactions also rejects single interaction events with 
extra (fake) vertices resulting from track reconstruction ambiguities. By comparing the fraction of events surviving 
the vertex cut with the fraction of single interaction events expected from the BBC cross section and the instantaneous 
luminosity, the vertex cut efficiency (fraction of single interaction events retained) is found to be 0.87 ± 0.02(syst) 
[0.90 ± 0.02(syst)]. In determining the DD gap fraction this efficiency is applied only to the total number of events, 
since the gap events have low central multiplicities and therefore are not likely to have fake vertices. 

The measured DD gap fractions, which are based on our experimental gap definition, Ar^^^ = rimax — Vmin, depend 
on the particle Et thresholds used. The correction factors needed to transform the measured gap fractions to gap 
fractions corresponding to the gap definition on which Eq. a is based, namely Ary" = In . 5/° ^ (In ^^ < 0, i — 1,2), 

were evaluated using the DD Monte Carlo simulation and found to be 0.81 [0.75] for y/s = 1800 [630] GeV. Correcting 
the measured DD gap fractions by these factors and for the vertex cut efficiency, and normalizing the results to our 
previously measured cross sections of 51.2 ± 1.7 mb [39.9 ± 1.2 mb] for events triggering the BBC arrays, we obtain 
2.51 ± O.Ol(stat) ± 0.08(norm)±0.58(bg) mb [2.16 ± O.Ol(stat) ± 0.06(norm)±0.65(bg) mb] for the DD cross section 
in the region A?/'^ > 3. 

The trigger acceptance, evaluated from the DD MC simulation, is 0.57±0.07(syst) [0.63±0.07)(syst)]. The uncer- 
tainty was estimated by considering variations in the simulation of small mass diffiaction dissociation. The acceptance 
corrected DD cross sections for A?7° > 3 are 4.43 ± 0.02(stat)±1.18(syst) mb [3.42 ± 0.01(stat)±1.09(syst) mb]. 

The corresponding cross sections predicted by Eq. 0, determined by the DD MC simulation, are 49. 4± 11.1 (syst) mb 
[27.7±5.5 (syst) mb], where the uncertainty is due to an assigned 10% systematic error in the triple-Pomeron coupling 
g(0) = K/3(0) [^. The ratio (discrepancy factor) of measured to predicted cross sections is Djjjj = 0.09 ± 0.03 
[0.12 ± 0.03], where the errors include the statistical and all systematic uncertainties. Recalling that Eq. correctly 
describes the shape of the Ai]'^ distribution, the deviation of D from unity represents a breakdown of factorization 
affecting only the overall normalization. This result is similar to that observed in SD 0,^2] , where the corresponding 
discrepancy factors, calculated from the fit parameters in Ref. [g|, are Dsd = 0.11 ± 0.01 [0.17 ± 0.02]. 

Our data are compared with the UA5 results Q in Fig. ??. The comparison is made for cross sections integrated 
over t and over all gaps of A77 > 3, corresponding to ^ = e^'^^ — 0.05 in SD. The extrapolation of our data from 
A77° > 3 (gaps overlapping 77 = 0) to Ar] > 3 (all gaps) was made using Eq. H and amounts to multiplying the Arj^ > 3 
cross sections by a factor of 1.43 (1.34) at y/s = 1800 [630] GeV, yielding 

aDoiVs = 1800 [630] GeV, A?7 > 3) = (4) 

6.32 ± 0.03(stat)±1.7(syst) mb 
[4.58 ± 0.02(stat)±1.5(syst) mb] 

The reported UA5 cross section values were obtained by extrapolating cross sections measured over limited large-gap 
regions to Arj > 3 using a Monte Carlo simulation in which the p and p dissociated independently with a (1/M^)e^* 
distribution llq] . For a meaningful comparison, we corrected the reported UA5 values by backtracking to the measured 
limited A77 regions using a (1/M'^)e^* dependence and then extrapolating to Ar] > 3 using Eq. 0. This correction 
increases the cross sections by a factor of 1.43 [1.19] at y/s — 200 [900] GeV. The solid curve in Fig. ?? was calculated 
using Eq. 0. The disagreement between this curve and the data represents the breakdown of factorization discussed 
above. The dashed curve represents the prediction of the renormalized gap probability model [|l2,n0|, in which the 
integral of the gap probability (first factor in Eq. ^ over all available phase space) is normalized to unity. The error 
bands around the curves are due to the 10% uncertainty in the triple-Pomeron coupling [y. Within the quoted 
uncertainties, the data are in agreement with the renormalized gap model. 

In conclusion, we have measured double diffiaction differential cross sections in pp collisions at -y/s =1800 and 
630 GeV and compared our results with data at •\/i=200 and 900 GeV and with predictions based on Regge theory 
and factorization. We find a factorization breakdown similar in magnitude to that observed in single diffiaction 
dissociation. The data are in agreement with the renormalized gap probability model [HOJ. 
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FIG. 1. Schematic diagram and event topology of a double diffractive interaction, in which a Pomeron (JP) is exchanged 
in a pp collision at center-of-mass energy ^/s producing diffractive masses Mi and M2 separated by a rapidity gap of width 
A?7 = i]max — rjmin- The shaded areas represent regions of particle production (mass and energy units are in GeV). 
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FIG. 2. The number of events as a function of rjmax and —r]min, the r) of the track or hit tower closest to »; = in 
the (anti)proton direction at ^/s = 1800 GeV: (a) data; (b, c, d) MC generated non-diffractive (ND), single- (SD) and 
double-diffractive (DD) events. The MC distributions are normalized by a fit to the data described in the text. 
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FIG. 3. The number of events as a function of l^rflxp ~ Vmax — ilmtn for data at ^/s = 1800 GeV (points), for double 
diffractive (DD) plus non-DD (MC) generated events (solid line), and for only non-DD MC events (dashed line). The non-DD 
events are a mixture of 97.3% non-diffractive and 2.7% single diffractive. 
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FIG. 4. The total double diffractive cross section for p{p) + p ^ Xi + X2 versus ^/s compared with predictions from Regge 
theory based on the triple-Pomeron amplitude and factorization (solid curve) and from the renormalized gap probability model 
(dashed curve). 
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